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Abstract- In this paper, we suggest a variant of Newton - Raphson method with which one can attain up to 
fourth order convergence provided the difference between f’(x) and /'(a) is negligible up to the desired 
accuracy, where 'a ' is the real root of /(a - ) = 0 . 
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I. INTRODUCTION 



The most widely used algorithm for solving a nonlinear equation in one variable i.e; 

f O)=0 (1.1) 

Where / \D czR -A R is a scalar function on an open interval D, by the use of value of function and its derivative 
is the well known Newton’s method given by 



f_0O 



(«= 0 , 1,2 ) 



(1.2) 



which has a quadratic convergence. 

As is known that derivation of Newton’s method involves an indefinite integral of the derivative of the function, and 
the relevant area is approximated by a rectangle. But, by approximating this indefinite integral by a trapezoid 
instead of a rectangle, Weerakoon and Fernando [4] suggested a Variant of Newton’s method which requires three 
functional evaluations at each step, given by 



= A . 

35+1 " [f’{x n ) + f’{x n+l )\ 

(1.3) 



(«= 0 , 1 , 2 ) 



where x =x 

n + 1 n 



n* n ) 



(1.4) 



which has a cubic rate of convergence. 

In this paper, we suggest a Variant of Newton-Raphson method which requires only two functional evaluations 
at each iteration apart from one extra evaluation initially. In section 2, we suggest the Variant of Newton-Raphson 
method whereas the rate of convergence of this method is obtained in section 3. In the concluding section, we give 
some numerical examples to exhibit the superiority of this Newton-Raphson’s variant. 
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II. VARIANT OF NEWTON - RAPHSON METHOD 



Let x Q be an initial approximation to the real root 'a' of the equation (1.1) is an open interval in which 
f(x) is continuous and let 

X* 0 =X 0 +h (2.1) 

and x l be the next two successive approximations. 

If x Q be in the vicinity of 'a' and if ‘h’ is small enough, then as derived in the Newton’s method, one can have 

/ K) 



*o=*o- 



and then the next approximate 'x l ' is to be obtained as 



x=x„- 



no 

A 

no 



1 "0 r U 



no 

by fixing f'( x 0 ) an d x 0 * s as given in (2.2). 

Again, keeping f\x ) unchanged, the next approximation x 2 to be obtained as 

no 



( 2 . 2 ) 



(2.3) 



x z= x x~ 



no 



(n= 0,1,2 ), 



(2.4) 



where x*=x— — — i— - . 

1 1 no 

In general, this variant of Newton - Raphson method can be defined as 

no 



x . =x - 

n+ 1 n 



no 



(«= 0 , 1,2 ) 



(2.5) 



no 



where X n —X n — — , which requires only two functional evaluations at each step leaving the one f'(x a ) 

f (O 

calculated initially. 

Algorithm 2.1: For a given x (} , compute the approximate solution x ] by the iterative scheme: 

* no , * n x B ) 

x ,—x where x —x 

71+1 " no " n no 

This algorithm requires only two functional evaluations at each step leaving the one derivative f'(x 0 ) 
calculated initially. Also, the efficiency index of the method \l~4 i.e the order of the method is four provided the 
difference between f'(a ) and /'(jt ) where 'a' the exact root of (1.1) be, is zero up to the desired accuracy 
neglecting evaluation of only once initially. 



III. RATE OF CONVERGENCE OF VARIANT OF NEWTON - RAPHSON METHOD 



Theorem 3.1 : Let a e D be a simple zero of the function f : D <^R — > R in an open interval D. If x Q is in the 

f ,^ 

verv vicinitv of the root ’a' then the algorithm (2.1) has fourth order convergence provided is nearly unitv 

no 

up to the desired accuracy. 

Proof: If 'a' be the exact solution of the eqn. (1.1), then 

/(«)= 0 (3.1) 

Let e n+1 and e be the errors at (n+ l/ h and n' stages and let X n+1 and x n be the («+ l) th and n' approximations to 
the root 'a ' of the eqn. (1.1). Therefore, we have 
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x . -e . +a 

n + 1 n +1 

x -e + a 

n n 



Now, f(x n ) = f(a + e n ) 



e „ 1 f"(a ) 2 1 /'"(a) 3 1 / A '(a) 4 sm 

— f (0 C > )\g H 6 H 6 H 6 + 0 (g 

" 2! f\a) n 3 \ f\a) n 4! /'(a) n 



f\a) a2 , /"'(«)„ 3 , ^' V («) e 4 +ofe 5 

72 ^72 



-e + - 



= f(cc) + f'(a)G + 

n 2! " 3! " 4! 

= ^'(a) K + + c 3 e* + c 4 e n 4 + 0(e/)] 



Where c = (/=2, 3,4...) 

7 /! / '(a) ' 

/(a- ) 

Now again, X n —X n n 



= a + e 



n* o) 
/•'(«) 
/■'Uo) 



[e + c.e 2 + ce 3 +c e 4 +0(e 5 )] (by 3.3 and 3.4) 

L 22 272 3 72 4 fl v 72 ' J 



= a + (1 - k)e - c„ke 2 -c.ke 3 - c, ./re 4 + 0(e 5 ), 

v ' n 2 n 3 n 4 n K n ' ’ 



where 



k — 



[M 

fW 



(3.2) 

(3.3) 



(3.4) 



(3.5) 

(3.6) 



f (x* n ) = f (a) + f’(a)[(l- k)e n - c 2 ke 2 -c 3 ke 3 ~c 4 ke n 4 +0(e s n )] 

+ ^^-[(1 - k) 2 e n 2 + c^ 2 e n 4 + 2{-Ar(l ~k)c 2 e 3 -k{ 1 - k)c 3 e n 4 } + 0{e s n )] 



2! 
/•"(«) 
3! 



[(1 - kfe 3 + 3{-/r(l ~kfc 2 e n 4 -k 2 { 1 ~k)c 2 e n 4 + 0(e;j )}] 



z* /F 

+ — [(l-^ ) 4 e„ 4 + 0 (^)] + . 



= /'(«)(! - /r)e„ -[/'(« )c 2 /r - ^^-(1 - ^) 2 ]e„ 2 
-[/'(a)c 3 A- +/"(a)A-(l-A-)c 2 - /"'(a)(l - Ar) 3 ]e„ 3 

-[/'(a)c 4 ^r + f"(a)k(l - k)c 3 + f ^ k(l- k) 2 c 2 + f ^ c 2 k 2 

-^ i< i-* ) 4 K‘ + o (e = ) 

= f'(«)(l-t)e, -[f(a)c 2 * - 2 ^ 1(1 -/r) 2 ]e„ 2 

- [f\a)c 3 k + f"{a)k{ 1 - k)c 2 -/"'(«)( 1 - ^) 3 ]e„ 3 

- [/-'(a)c 4 /r + f\a)k{{l - k)c 3 + -|-)+ — ^*( 1 - £) 2 c 2 

_/^a2 ( i _^ )4] ^ 4+0(e 5 ) 

4 ! 

Now, 
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xl ~ f ^ n = a + (1 ~k)e n -c 2 ke n 2 -c 3 ke n 3 - c 4 ke n 4 






- ^i(i - + { ^c 2 k - rw-v - k f ]e 2 

f\x 0 y n f\x 0 ) 2 2 ffx 0 y 

+ [^ + ^(l_,) C 2 _m ( l_,) 3 ] e 3 
f\x 0 ) 3 /'(*<,) f, (x 0 ) 

/'(a) „ f , , /"(«) , c\k /'"(a) ,, 2 . 



+ [ — - — —c.k + — - — —A - ( (1 - k )c, + -*— ) + 

L f'f „ ^ 4 /* O' A VV y 3 o' 



n* 0 ) 



(l-^) 4 ]e n 4 + 0 ( e /) 



2 /'On) 



A(1 - A) z c 2 



By (2.5), we have 

X n + 1 = « + C 1 “ *) e fl - C 2^« 2 - C 3 ke n ~ C <t ke n 

- *( 1 - ■ *>, + [<- 2 * 2 - iTyffjC 1 " fr ) 2 ]^, 2 

+ 1 \CJ 2 + 7^*(1 - * )c 2 - ^ 2 (1 _ k f ]e / 

4 4+ s 

4 !/'(* 0 ) " " 

Where 'A -1 is as given in (3.6). 

= « + (1 - k ?e n ~ [*(1 - A)c 2 + /" (a) (1 - A) 2 ]e 2 

n ,, , . /"(a) , . , . /'"(a) . , . 3 , 3 

- [A(l - A)c, — — A(1 - A)c 9 + — ^-(1 - A rle„ 

3 a'(^o) 2 n * Q ) 



[A(l - A)c. 



/•"(«) 



A((l - A)Co H 7 - ) 



^ 2 * /■"'(«) 



2 /'K) 



A(1 - kyc 7 



+ / (Qf) (1- A) 4 ]e 4 +6>(e„ 5 ) 

4 !/'(* 0 ) " " 

Therefore, by (3.8), (3.2) the method (2.5) takes the form 
a + e « + i = a + ( 1 “ Ar ) 2e n ~[A( 1 ~A)c 2 + ^r^y (l -A) 2 ]e fl 2 

- [A(l - A)c g - f (a) A(1 -A)c 2 + 1 (ar) (1 - A) 3 ]e 3 
3 f\x 0 ) 2 f(x 0 y 

f\ot) ^ /"(«) _ 2 . 



A(1 - A) c 2 



^+1 = (1 - KYy - [d - + TTTT^rd - kf ]e n 2 



- [A(l - A)c, ^2_A((1 - k)c, + -^— ) - 

L ^ ^ 4 /'(a- 0 ) 3 2 2 2f\x 0 ) 

+ f ‘ V (a) ( 1 - A) 4 ]e 4 + 0 (e„ 5 ) 

4 !/'(^ 0 ) 

2 „ m /, \„ /, , ^ ) 

2f\x 0 ) ' " 

-[(1 -k)c 3 k- f (a) A(l-A)c 2 + - (a) (1 - k ) 3 ]e 3 
-[(1- A)c 4 A - / ^ A((l - A)c, + C -l—)~ f ( - a) k{l-k) 2 c 2 

fiVf \ n* 0 ) 3 2 2/'(* 0 ) 

+ f (a) (1 — /r ) 4 ]e 4 +0(e n 5 ) 

4!f(x a ) " 71 



(3.8) 



(3.9) 



(3.10) 
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It can be seen from the above step that the terms up to and including e 2 vanish if ‘ k ’ is nearly unity up to the desired 
accuracy and if it does so, we have 



e 



n+l 



f\a) 

2 f\* o) 



e>0(e 5 J 



(3.11) 



That means e n+1 cc (3.12) 

The method (2.5) has a fourth order convergence if k=\ . 



IV. NUMERICAL EXAMPLES 

We consider few numerical examples considered by Weerakoon and Fernando [4] and by Grewal [1] and the 
method (2.5) is compared with the methods (1.2) and (1.3). The computational results are tabulated below and the 

results are correct up to an error less than 5 x 1CT 7 . 

Table 1 



Function 


*0 


i 


NOFE 


Root 


f{x) 




NM VNM VNR 


NM VNM VNR 




(1) £?* -1 +x-3 


1.4 

1.5 


3 2 2 

3 2 2 


6 6 5 

6 6 5 


1.442854 


(2) 2*-log 10 ;r-7 


3 

3.1 


3 2 2 

3 2 2 


6 6 5 

6 6 5 


1.404492 


(3) x log 10 x -1.2 


2.5 

2.6 


3 2 2 

3 2 2 


6 6 5 

6 6 5 


0.25753028 


(4) cos(a')-* 


0.7 

0.8 


2 2 2 

3 2 2 


4 6 5 

6 6 5 


0.7390851 


(5) x 2 -e x -3x + 2 


0.3 

0.4 


2 2 2 

3 2 2 


4 6 5 

6 6 5 


0.2575303 



NM - Newton’s Method VNM - Variant of Newton’s Method 

VNR - Variant of Newton - Raphson NOFE - Number of Function Evaluations 

i-Number of iterations to approximate the root to 7 decimal places 



V. CONCLUSION 

From the above tabulated values, we conclude that the method (2.5) has an advantage of more rapid convergence 
with only two functional evaluations from second iteration onwards, over the other methods discussed in this paper. 
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